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We perform a detailed analysis of how high-order harmonic generation (HHG) and above-threshold 
ionization (ATI) can be controlled by a time-delayed attosecond-pulse train superposed to a strong, 
near-infrared laser field. In particular we show that the high-harmonic and photoelectron intensities, 
the high-harmonic plateau structure and cutoff energies, and the ATI angular distributions can be 
manipulated by changing this delay. This is a direct consequence of the fact that the attosecond pulse 
train can be employed as a tool for constraining the instant an electronic wave packet is ejected in 
the continuum. A change in such initial conditions strongly affects its subsequent motion in the laser 
field, and thus HHG and ATI. In our studies, we employ the Strong-Field Approximation and explain 
the features observed in terms of interference effects between various electron quantum orbits. Our 
results are in agreement with recent experimental findings and theoretical studies employing purely 
numerical methods. 
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I. INTRODUCTION 

High-frequency light pulses of attosecond (10~^^s) du- 
ration have caused a breakthrough in metrology, allowing 
one to resolve and control dynamical processes occurring 
in the subfemtosecond and subangstrom scale [3| . This is 
in particular possible due to the fact that one attosecond 
is roughly the time it takes for light to travel through 
atomic distances. Attosecond pulses may be used, for 
instance, for tracing the motion of bound electrons 0, 
exciting inner shell electrons la] , controlling molecular 
motion Ij, electron emission o\, high- harmonic genera- 
tion ^ Jj or above-threshold ionization 8]. 

Such pulses owe their existence to the fact that high- 
order harmonics, which are generated by rare gases in- 
teracting with intense (/ ~ lO^^W/cm^), near- infrared 
(IR) laser fields of femtosecond duration, are nearly phase 
locked. Hence, by superposing a large number of har- 
monics one may obtain attosecond pulses in the extreme 
ultraviolet range There exist two approaches to 

attosecond-pulse production. One may generate isolated 
attosecond pulses from few-cycle laser pulses |^ 9, 10], 
or attosecond- puls e trains from laser pulses comprising 
several cycles [ll|. Such attosecond pulse trains have 
been predicted in the late 1990's and experimen- 
tally realized a few years later . Since then, there has 
been considerable improvement in their contrast, dura- 
tion and peak intensity. In particular, an intrinsic chirp 
of the attosecond emission was evidenced, setting a lower 
limit on the duration of the attosecond pulses obtained by 
grouping an increasing number of harmonic orders |l3j] . 
Thereby, important issues are an adequate choice of the 
group of harmonics, and optimal propagation conditions 
in the gaseous media. In this context, it is of interest to 
control the atto-chirp and more generally the intensity- 



and order-dependent phase of the atomic dipole 0] . 

Recently, it has been shown that an attosecond-pulse 
train superposed to an intense laser field can be used 
for manipulating high-order harmonic generation (HHG) 
0, 13 and above-threshold ionization (ATI) Q . In par- 
ticular, the resolution, intensity and maximal energies of 
both ATI and HHG spectra turned out to be strongly 
dependent on the time delay between the laser field and 
the attosecond pulse train. This behavior has been pre- 
dicted theoretically ^ ^] and observed experimentally in 
both cases 0, ■ It has been interpreted in the light of 
the physical mechanisms governing both phenomena. 

In HHG, an electron reaches the continuum by tun- 
neling or multiphoton ionization at a time t' , propagates 
in the continuum and, subsequently, at a time t, recom- 
bines with its parent ion 15]. In this process, the kinetic 
energy acquired by this electron is released in form of 
high-frequency radiation. Similarly, for ATI, an electron 
leaves an atom, either reaching the detector, or elastically 
recolliding with it. This leads to the low- or high-order 
ATI peaks, respectively. The attosecond pulse train can 
be used for controlling the ejection of such an electron 
in the continuum. Specifically, the time t' can be fixed 
by changing the time delay between the attosecond pulse 
train and the IR field. This constraint upon the ioniza- 
tion time can be used as a tool for selecting a particular 
orbit of an electron reaching the continuum JJ , or a 
particular momentum transfer from the IR field to the 
electron . As a direct consequence, the ATI and HHG 
spectra can be manipulated. 

In this paper, we investigate the effects reported in 
0, Q for HHG and ATI in detail, using a quantum- 
mechanical formulation of the above-stated processes 
|l6l [T7|. | l 8ll. within the Strong-Field Approximation 
(SFA) lis]. This approach has some advantages over 



2 



that employed in I3, i.e., the numerical solution of 
the time-dependent Schrodinger equation (TDSE). In- 
deed, even though the SFA involves a higher degree of 
modeling and physical approximations, it is much less 
demanding, from the computational viewpoint. This is 
particularly important for a realistic description of the 
macroscopic response of the gaseous system, in which 
the single-atom response serves as input for propagation 
codes. Furthermore, it provides a transparent and clear 
physical interpretation, in terms of the orbits of an elec- 
tron leaving or recombining with its parent ion. This is 
a further advantage over the TDSE, for which it may be 
difficult to extract the physical mechanisms involved. 

In this work, we put a particularly strong emphasis on 
the quantum interference between the different possible 
paths of an electron ejected by the attosecond-pulse train. 
Specifically, we investigate how the initial conditions with 
which the electron reaches the continuum influence such 
an interference, and which consequences this has on the 
ATI or HHG spectra. We show that the attosecond pulses 
can be used for controllingthis interference, and explain 
the features observed in Q, ll| in terms of electron 
orbits. 

This paper is organized as follows: In SecE] we briefly 
discuss the SFA for ATI and HHG, with emphasis on the 
changes introduced by the attosecond-pulse train. Subse- 
quently (Sec. IIII|I . we show how several features of both 
phenomena are affected by it, and provide physical ex- 
planations. Finally, in Sec. IIVI we give a summary of 
such results and state the main conclusions of this work. 



II. TRANSITION AMPLITUDES 

A. Above-threshold ionization 

For ATI, we will concentrate on the direct electrons, 
which reach the detector without rcscattering with their 
parent ion. In this case, the ATI transition amplitude 
reads 

bp{t) - z /* E{t')dt'd,{p + A(t'))e^^(P'*'*'\ (1) 

JQ 

with the action 



Sip, 



(2) 



Eq. can be understood as follows: At an instant 
t', an electron is ejected by the external laser field from 
the atomic state in which it is initially bound, reaching a 
continuum state. It then remains in the continuum until 
a time t, when it reaches the detector with final momen- 
tum p. In Eq. lUJ, A{t), Ip, and E(i') = -dA{t')/dt' 
denote the vector potential, the atomic ionization poten- 
tial, and the electric field at the time of its release, respec- 
tively. We take the dipole matrix element dz{p + A{t)) = 



(p -I- A{t) \ Y.€z IV'o) for an initial Is state, where the vec- 
tor denotes the polarization axis. Specifically, this 
yields 



2i9/4r 



5/4 ' 



pe2 



[p2 + 2Jp 



(3) 



We will now write the transition amplitudes for the 
specific case in which an attosecond-pulse train is present. 
We assume that this train is composed of an infinite num- 
ber of odd harmonics. This yields 



Eatto(0 — EhTT 



n=0 



Mr 

a{t) 



o(t 



(4) 



where w, Eh, and cr(i) denote the laser field frequency, 
the attosecond-pulse strength and the train temporal en- 
velope, respectively. In this paper, we assume that the 
envelope is either given by 



(T(t) -exp[C|t|], 



(5) 



where C ^ 1, or by cr = const. Physically, the latter case 
corresponds to an infinitely long attosecond-pulse train. 
We approximate the laser field by a monochromatic wave 
E;(t) = —dAi{t)/dt, with the vector potential A/(t) — 
2 ,yUp cos{ujt — (j))ez of frequency uo dephased of (f) with 
respect to the attosecond-pulse train. 

We now assume that the attosecond pulse train re- 
leases the electron in the continuum and that its subse- 
quent propagation is determined only by the monochro- 
matic field. This implies that E(t') ~ Eatto(^') and 
A(t) ~ Ai{t) in Eq. For the direct ATI electrons, 
such assumptions lead to the transition amplitude 

bp ^ [d cos a — id^ sa\a\x{cS°\ (,11 / uj) 

-f [d_ cos a — id+ sin Q;]x(a*-'^\ C"'/'^); (6) 

where (i± = d2(p-|-A(t'))±(i2(p— A(i')), 6* gives the angle 
between the laser field and the final electron momentum, 
and the argument a = 2p ^JTJ^ / lo cos 9 sn\ 4) comes from 
the non-trivial part of the action. Furthermore, in Eq. 



x(a,C7r/w) 



sinh(C7r/w) 



■ cos(a) -I- cosh(C7r/cj) 



(7) 



a(°' = [Ip + Up+ - uj)t:/oj and a^^^ = [Ip + Up + 
/2)it/uj, respectively. Eq. O exhibits sharp peaks at 

cos a = 1 (i.e., for a = 2mT). Specifically, a^"' and a^'^^ 

lead to expressions with steep maxima around the odd 

and even ATI peaks, respectively. 

In the limit of infinitely long attosecond-pulse trains. 



6p°'' , where 



K""^ ^ X! ~ id+sina] 6{^ + Up+Ip-2nuj) (8) 
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and 



COS a 



id. 



. sin a] S{ Y+Up+Ip-{2n+l)u}) 



(9) 

gives the contributions from the even and odd ATI peaks, 
respectively In Eq. JHJ and ©, one can explicitly 

notice that the electron leaves with vanishing momen- 
tum if 2nu! = {Ip + Up). Physically, this means that the 
ionization potential Ip is effectively shifted by the pon- 
deromotive potential Up. Such an effect is noticeable for 
relatively high driving- field intensities. 

In some of the subsequent results, in order to perform 
a more direct comparison with the results existing in the 
literature, we will integrate the photoelectron yield over 
the solid angle. In this case, the yield will be given by 



r]{p) = 27ry \bp\'^d{cosd), 



(10) 



where bp is given by Eq. ®. 

Furthermore, apart from initial Is states, we will con- 
sider the so-called broad gaussian (GBR) approximation, 
for which ^^(p -|- A{t)) = const. Physically, this ap- 
proximation means that the electron is initially in an ex- 
tremely localized bound state (for details see, e.g., [l^'). 
In this particular case, and in the limit of infinitely long 
attosecond-pulse trains, the integral ((Tn|) can be solved 
analytically and gives 



vip) 



1 ± 



sin(2Q;o) 
2ao 



(11) 



where ao — 2pyJU^/u} suvcj)^ and the positive and neg- 
ative signs correspond to the odd and even ATI peaks, 
respectively. 



B. High- harmonic generation 



We will, however, adopt a slightly different convention: 
instead of considering the time delay in the infra-red field, 
we will take A;(t) = 2yJU^cos{u)t)ez and insert the delay 
id — 'i'l^ in the attosecond pulses. Such a convention will 
facilitate the discussions in Sec. III.B, in terms of electron 
orbits. Under such assumptions, Eq. (|12|l reads 



hn = 



— Yi-ir / dt / d3ytexp[z5(t,i'„,r!,k)] 
^ ^0 J-^J 



d:{-k + Ai{t))d,{k + Ai{t'J). 



(14) 



In this case, the ionization time is being fixed at t'j = 
td + mr/uj, with n an integer, because the electron is 
being released into the continuum at such an instant by 
the attosecond pulses. 

The transition amplitude is further simplified in the 
sense that only the integrals in k and t must be solved. 
We evaluate this equation employing saddle-point meth- 
ods. For that purpose, we consider k and t so that 
S{t,t',n,'k) is stationary, i.e., that dS/dls. = and 
dS/dt = 0. This yields the saddle-point equations 



and 



dT [k + Ai(T)] =0, 



2{n -Ip) = [h + Ai{t)]' 



(15) 



(16) 



whose solutions can be directly associated to the classical 
orbits of an electron returning to [Eg. H15() ] and recom- 
bining with [Eq. (|16|l ] its parent ion 18] . In the presence 
of the attosecond-pulse train the saddle-point equations 
for t and k can be combined, so that the electron return 
time is given by 



The HHG transition amplitude is given by an expres- 
sion which is slightly different from Eq. namely 

bn =i dt dt' / d^kdl{k + A{t))d,{k + A{t')) 

J — ooJ — oo J 

E(t')exp[iS{t,t',n,k)], (12) 
with the action 

1 r* 

s{t,t',n,k) ^ -- / [k+A{T)fdT-ip{t-t')+nt. (13) 

2 Jt' 

Eq. El describes a physical process in which an electron 
is freed at a time t' , propagates in the continuum with 
momentum k from t' to t, and, at this time, recombines 
with its parent ion, generating a harmonic of frequency 

n. 

We will now apply the same physical assumptions as 
in the previous sections to the transition amplitude (|12|l . 



sinujt - (-1)" sinLjtd (17) 
= [Lu{t- td) - mr)] (cos tut T ^ {n - Ip) / {2Up)^ . 

For details on the approximation employed we refer to 
Ref. 113. 

Alternatively, one may follow the procedure in 
Ref. ^3 , and compute the Fast Fourier transform for the 
time-dependent dipole moment in the Strong-Field Ap- 
proximation. Such an approach involves a single saddle 
point approximation in the intermediate electron momen- 
tum, and two numerical integrations. Hence, it possesses 
the advantage of taking into account all sets of orbits, 
and not only the above-mentioned pairs. The procedure 
adopted in this paper, however, allows a more detailed 
assessment of quantum-interference effects. For compar- 
ison, we provide the time-dependent equations employed 
in JM ' together with a brief discussion, in the appendix. 
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III. RESULTS 
A. Above-threshold ionization 

As a starting point, we will investigate the effect re- 
ported in 8j, in which the yield extends to much higher 
energies for vanishing delay. In @, this is justified by the 
following argument: if the electron leaves the atom at the 
zero of the IR field (i.e., at = mr), its vector potential 
is at is maximum, i.e., A(t') = -^I^JTTp. Hence, when the 
electron reaches the continuum, under these conditions it 
acquires maximum possible drift momentum. Therefore, 
it reaches the detector with higher energies. Since this 
effect only influences how the electron is ejected in the 
continuum, this should affect both direct and rescattered 
electrons equally. Therefore, we expect that Eq. © and 
©, for the direct electrons, exhibit such a behavior. 

Fig- n displays ATI spectra obtained from Eq. ©, for 
time delays < < 27r. In order to perform a direct 
comparison with the results in @ , we integrate the yield 
over the solid angle [Eq. H1U|| ] . In the upper part of the 
figure [Fig. n(a)], we consider the same parameters as 
in [sj, while in the lower part [Fig . n](b)] we choose a 
higher intensity for the IR field |2^ . In the figure, there 
exist two distinct behaviors, within two energy regions. 
For photoelectron energies near the low-energy end of 
the spectra, the peaks exhibit maximal intensities near 
= (2n-|-l)7r/2, i.e., if the electron is released close to the 
times for which the electric field is maximal. For higher 
energies, the yield presents a clear minimum around such 
a delay, and the spectra extends to much higher energy 
ranges if the delay is vanishing, i.e., if the electron reaches 
the continuum at the electric-field crossings. This be- 
havior agrees with the findings in , and becomes more 
extreme as the laser-field intensity increases. 

Both energy regions appear to be dependent on the 
driving- field intensity. For instance, for Fig.n(a), max- 
ima near = (2n + l)i'/2 are observed roughly for 
ATI orders 12 < iV < 20 (0 < pV2 < IS.SeF), while 
for Fig. □(b), this region lies within 32 < < 40 
(0 < pV2 < 16.4ey). Beyond such energy regions, the 
yield exhibit minima for such delays. Apart from that, for 
high-energy photoelectrons, there exist additional min- 
ima for the higher driving- field intensity [c.f. Fig.n](b)]. 

Such a behavior can be readily seen from Eqs. Q and 
(|^. or, more specifically, from the angle- integrated yield 
[Ea. l(Tm) ]. Indeed, we expect to be the superposition 
of two types of contributions, coming respectively from 
the non-trivial part of the action a and from the dipole 
matrix elements (i±. Physically, the former corresponds 
to the phase the electron acquired leaving at a particular 
instant, and the latter are specific to the initial state in 
which the electron is bound. 

In Fig. 13 we display the isolated contributions from 
such terms, integrated over the solid angle [Eq. H10|) ]. for 
distinct ATI orders. In order to disentangle the effect due 
to the dipole matrix element and that due to the non- 
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FIG. 1: (Color online) Direct above-threshold ionization 
spectra, integrated over the solid angle and in linear scale 
(Eq. not '), for argon (/p = 0.58 a.u.) interacting with a 
monochromatic field of frequency lj = 0.057 a.u. and an 
attosecond-pulse train (Eq. (|U), as functions of the time 
delay between both driving fields. Panels (a) and (b) have 
been computed for an intensity of J = 3.3 x 10^"^W/cm^ and 
J = 5 X lO^^W/cm^ for the low-frequency laser field, respec- 
tively. In panels (a) and (b), the ATI orders (labelled with 
the number of photons from the ground state) correspond to 
photoelectron energies roughly between and 36.5 eV, and 
and 47.4 eV, respectively, due to the ponderomotive shift. 
The yield has been computed from Eq. with C, = 0.04. 
The intensities of the attosecond pulse trains have been cho- 
sen as //i = //lO, i.e., one tenth of those of the low-frequency 
field. 



trivial part of the action, in addition to considering Is 
initial states, we computed spectra in the broad gaussian 
(GBR) approximation, for which (iz(p -I- A(t)) = const 
[Eq. Hll|) ]. In this case, d+ = 1 and d- = 0, so that the 
effect caused by the non-trivial part of the action, i.e., 
sin a and cos a in Eq. (6), can be directly visualized (c.f. 
gray curves in the upper panels). 

We observe that the non-trivial part of the action leads 
to a highly oscillating pattern, which, for even or odd ATI 
peaks, have maxima or minima at = titt, respectively. 
As the electron energy increases, the oscillating region 
becomes increasingly flatter, and the amplitude of the 
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FIG. 2: Upper panels: modulus squared of the transition am- 
plitudes 10 and integrated over the solid angle, for Is ini- 
tial states (black line) and within the broad gaussian approx- 
imation (grey line), as functions of the time delay between 
the attosecond-pulse train and the IR field. Lower panels: 
|d+P and integrated over the solid angle, for Is initial 

states. We considered argon under the influence of an IR fleld 
of / = 3.3 X lO^^W/cm^ superposed to an attosecond-pulse 
train of /h = J/10 (Panel (a) in Fig. 2). The curves in the 
upper panels have been normalized so that their maximum 
values are unity, and those in the lower panel to the maxi- 
mum value of D\~\ The numbers N in the figure denote the 
ATI orders. 



oscillations are always larger near (f) — mr. It is, however, 
clear that the action, by itself, is not responsible for the 
behavior exhibited in Fig.nand shown in detail in Fig. 12 
for initial Is states (c.f. black curves in the upper panels). 
In fact, even though this term leads to several extrema, 
it does not cause an overall enhancement or suppression 
of the yield near = titt, as observed for the low- and 
high-order ATI peaks respectively. 
The integral 

D^f"^ ^2tt J \d±\'^d{cos0) (18) 

of over the solid angle sheds some light in this be- 
havior (c.f. lower panels in Fig. |2J). The contributions 
from d_ exhibit a minimum at (/) = (2n+l)7r/2, while the 
behavior of those from (i+ depend on the energy region in 
question. Near the ionization threshold, the latter contri- 
butions exhibit maxima near the delays (f> = (2n-f l)7r/2. 
As the ATI order increases, such maxima become flat- 
ter, until the concavity of Dj'^'^ changes and the yield 
starts to exhibit maxima dA, (f) — mr. Such contribu- 
tions are dominant, with respect to those from d-, and 
are present both in Eq. JHJ and 0. This behavior gets 
more pronounced as the driving-field intensity increases 
(higher intensities lead also to sets of secondary maxima 



FIG. 3: (Color online) Angular distributions for the ATI di- 
rect electrons, in a logarithmic scale, for the same field and 
atomic parameters as in Fig. 1. The upper and lower panels 
have been computed for an intensity of / = 3.3 x 10^'^W/cm^ 
and J = 5 X 10^*W/cm^ for the low-frequency laser field, re- 
spectively. The yields are displayed as functions of the time 
delay </!> — ujtd and orientation angle 6 between the laser field 
and the final electron momenta. The ATI orders N are indi- 
cated on the upper parts of the panels, and all contours have 
been normalized to the highest yield in the figure (left panels). 



in Dj , which contribute to the additional structure ob- 
served in Fig.n](b)). One should also note that, very near 
the threshold, for (f> = m:, d'j'^ > d\'^\ This causes the 
minima observed for ri{p) in Fig. 1(a), for initial Is states. 
In contrast, the yield from the broad-gaussian approxi- 
mation exhibits maxima in this case. 

Therefore, the behavior reported in Q is not only due 
to the phase the electron obtains when it is ejected in 
the continuum by the attosecond-pulse train, or, in other 
words, to the momentum transfer from the field to the 
electron. Additionally, it is related to the initial state in 
which the electron is bound, which, in our model, is a 
Is state. Indeed, such a behavior is absent in the broad 
gaussian approximation, which, physically, corresponds 
to an electron initially bound in a highly localized state. 
Such a strong dependence on the initial electronic bound 
state is, in a sense, quite unusual for phenomena occur- 
ring in the context of atoms in strong laser fields. 

The ATI angular distributions, depicted in Fig.|21 con- 
firm this trend. In fact, the figure shows distributions 
symmetric upon 9 n — 9 and (/) — > tt — </), in which it 
is possible to identify three distinct energy regions. Near 
the threshold {N =15 and iV = 32 for the lower and 
higher intensity, respectively), there exist annular max- 
ima, centered around (j) = 7r/2. Such maxima are most 
pronounced for 9 = and 9 — n, but also occur for other 
alignment angles. As the photoelectron energy increases 
(middle and right panels in the figure), the ring-shaped 
maxima split and move towards 4> = nir and 9 = mr, 
which are the maxima observed for high enough ATI or- 
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ders. Other, less pronounced maxima are also seen, as 
well as additional substructure. Such a pattern is due to 
interference effects and comes from the non-trivial part 
of the action. As the driving-field intensity increases, this 
pattern gets more complicated, but the overall behavior 
remains. 



B. High- harmonic generation 

Subsequently, we address the question of whether the 
effects reported in 0, , for HHG in a near-infrared laser 
field superposed to a time-delayed attosecond-pulse train, 
are present in our framework. Specifically, intensity vari- 
ations from one to two orders of magnitude in the plateau 
harmonics have been observed, whose prominence, and 
even existence, was strongly affected by the time delay 
(j) = Lutd- Moreover, the resolution of the harmonic peaks, 
and the plateau structure, with enhancements of either 
all plateau harmonics, or only near its low-energy end, 
were observed to depend very much on this parameter. 

In 0, such effects have been attributed to quantum- 
trajectory selection: In HHG, in the absence of the at- 
tosecond pulses, there exists one set of trajectories, which 
mainly contribute to the yield. Such trajectories merge 
at the cutoff, i.e., the maximal energy for which HHG oc- 
curs, and correspond to a longer and a shorter excursion 
time for the electron in the continuum [T^ . By using an 
attosecond-pulse train, one is splitting this set into two, 
due to the fact that the electron is released in the contin- 
uum with non-vanishing velocity, and, at the same time, 
selecting a trajectory of this set, by fixing the electron 
emission time. 

Thus, the attosecond pulses allow one to perform 
quantum-trajectory selection at the single-atom response 
level. For solely an IR laser field, it is only possible to per- 
form such a selection at the macroscopic level, by control- 
ling the propagation conditions p^ . In support of this 
argument, classical simulations have been performed in 
for an electron returning to its parent ion. Thereby, 
the initial electron velocity has been estimated from the 
average frequency in the attosecond-pulse train, which 
in has been composed by five harmonics (N = 11 to 
N = 19). 

In Fig. 21 we present HHG spectra computed with Eq. 
(|14|l . in the broad gaussian approximation. Such spec- 
tra are highly dependent on the time delay td between 
the low-frequency field and the attosecond-pulse train, 
and differ in the cutoff energy, and in the high-harmonic 
strengths and resolutions. Specifically, in the exam- 
ples provided in the figure, the cutoff energies vary from 
roughly Ip + l.SUp [Fig. |l(a)] to Ip + 2.Wp [Fig.|l(c)], 
and the intensities of the plateau harmonics differ in up 
to two orders of magnitude. Furthermore, for partic- 
ular time delays, the plateau harmonics are unequally 
enhanced. For instance, in Figs.0](b) and21(c), the low- 
plateau harmonics are more than one order of magnitude 
larger than those near its high-energy end. Such a group 
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FIG. 4: High-harmonic spectra for neon (Ip = 0.79 a.u.) in- 
teracting with a monochromatic field of intensity / = 5 x 
lO^^W/cm^ and frequency uj = 0.057 a.u., and a train of at- 
tosecond pulses of intensity Jh = lO^'^W/cm'^ (Eq. (QJ), and 
several time delays between the attosecond pulse train and 
the IR field. We considered the orbits with the five shortest 
pairs of excursion times t — t'. In the figure, we only show the 
harmonic orders for which Q > Ip. 



varies from the lower half [Fig. 2|(b)] to the lowest few 
harmonics [Fig. El(c)] of the plateau. The spectra are 
identical over ojtd — > ujtd ± tt (not shown). This is a con- 
sequence of the fact that \dfi{ujtd)f = |dn(wtd ±7r)|^. In 
particular, for Lutd — m:, the plateau harmonics exhibit 
maximal intensities in this context. 

The effects observed in Fig.J^are in qualitative agree- 
ment with those reported in [a^Q- In order to find out 
whether they are due to quantum-orbit selection, they 
will be analyzed in terms of the orbits of an electron re- 
combining with its parent ion. For that purpose, we solve 
Eq. (17) for a fixed start time. The real parts of its solu- 
tions can be directly associated to the return times t of 
a classical electron in an external field. 

Fig. [3 depicts the real parts of such return times for 
the six shortest orbits, as functions of the harmonic order, 
for the time delays in the previous figure. The orbits are 
denoted by the numbers which increase with the 

electron excursion time in the continuum. Depending on 
the time delay, the return times, and the maximal har- 
monic energies are significantly altered. This shows that 
the changes in the initial conditions of the electron being 
ejected considerably affect its subsequent dynamics. Fur- 
thermore, the return times occur in pairs, i.e., for a single 
ionization time, there is always a shorter and a longer 
travel time for the electron in the continuum. Such pairs 
of orbits coalesce at the maximally allowed harmonic en- 
ergies, which considerably vary with the delay (f). Such 
variations are explicitly given in Table I. In particular. 
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FIG. 5: (Color online) Real part of the recombination times 
as a function of the high-harmonic order, for neon (Jp = 0.79 
a.u.) interacting with a monochromatic field of intensity 
/ = 5 X 10^*W/cm^, and frequency uj = 0.057 a.u. and an 
attosecond-pulse train of intensity Jh — lO^^W/cm^, com- 
posed of an infinite number of harmonics (Eq. (HJ). Parts 
(a), (b), (c) and (d) correspond to delays (f) — Q, (j) = 0.257r, 
(j) = O.Stt and 4> = 0.757r between the attosecond-pulse train 
and the low-frequency driving wave, respectively. In the fig- 
ure, we only show the harmonic orders for which fl > Ip, and 
the pairs of orbits are indicated by the natural numbers (i, j), 
which increase according to the electron excursion time in the 
continuum. The ionization times are indicated by the blue 
dashed lines, and the numbers in the upper part of the fig- 
ure give the approximate kinetic energy of the electron upon 
return, in units of the ponderomotive energy Up. 



Orbits Ei^in{t,t')/Up 

(j> = = O.257r (I> = 0..5tt = 0.757r 

(1,2) 1.80 0.93 0.26 

(3,4) 1.65 2.11 2.50 2.43 

(5,6) 1.62 L37 L25 1.32 



TABLE I: Maximum kinetic energies Ey^in(t,t') — {Nuj — 
Ip)/Up, in units of the ponderomotive energy, for an elec- 
tron along the three pairs of orbits displayed in Fig. |^ No 
entry means that this energy is very close to the ionization 
potential. 



the largest variations occur for the set of orbits (1,2), 
namely from the ionization potential to Ip + l.SUp. 

Each of such local maxima potentially corresponds to a 
cutoff in the harmonic spectra. Whether, however, they 
will lead to only minor features or to a sharp decrease 
in the high-harmonic yield will depend on three main is- 
sues. First, it is necessary that the electron reaches the 
continuum with a large ionization probability. Second, 
due to the fact that the electronic wave packet spreads 
with time, it is desirable that such a wave packet spends 



as little time as possible in the continuum. This yields 
a large overlap between such a wave packet and that of 
a bound electron, and, consequently, prominent harmon- 
ics. Finally, there may be effects related to quantum in- 
terference, which suppress or enhance particular sets of 
harmonics. Since, due to the attosecond-pulse train, the 
electron is being ejected with a large probability for all 
cases, we expect that the main difference will be due to 
wave-packet spreading. Hence, the shorter the electron 
excursion time t = t — t' in the continuum is, the more 
prominent the corresponding harmonics should be. 

We will now analyze the contributions of each pair of 
orbits to the spectrum. Such contributions are presented 
in Fig. El Panel|Sl(a) shows that, for an attosecond-pulse 
train in phase with the low-frequency field, the main con- 
tributions come from the orbits (1,2). This occurs for two 
main reasons. First, the harmonics from this set of or- 
bits are two orders of magnitude larger than those from 
the remaining pairs. This is a consequence of the fact 
that T is very short for such a pair. Furthermore, for 
this specific phase, the cutoff energy for (1,2) extends to 
a relatively large harmonic order (c.f. Fig. [Sjand Table 
1), so that the contributions from the remaining pairs of 
orbits are overshadowed. 

For a delay of </> = 0.257r [Fig.El(b)], this picture starts 
to change. In fact, although the harmonics from (1,2) are 
much stronger than those from the remaining pairs, the 
maximal kinetic energy for (1, 2) lies near a much lower 
harmonic order (near N = 35). Thus, beyond such en- 
ergy, the corresponding contributions are exponentially 
decaying and the harmonics intensities, as well as the cut- 
off energy near 2.11Up, are mainly determined by (3,4). 
This causes the double plateau structure in the figure, 
and provides an alternative explanation for the enhance- 
ments only in the low-plateau harmonics reported in j^. 

For (j) — O.Stt [Fig. El (c)], instead, there is a very steep 
intensity drop at the the low-energy end of the plateau. 
This is a consequence of the very low cutoff for the orbits 
(1,2), near the 19**^ harmonic. For higher energies, the 
main behavior in the spectrum is determined by (3,4). 
Finally, for (j) = 0.757r [Fig. |Sl(d)], the maximal energy 
for (1, 2) is so close to the ionization potential, that such 
a pair does not play any role. Consequently, there is an 
overall decrease of at least one order of magnitude in the 
spectrum. At the same time, the orbits (3,4) become 
increasingly shorter. Around Lutd = 0.757r, such latter 
orbits have become the shortest contributing and the in- 
tensity starts to increase, until, at ujtd = tt, the same 
pattern as for vanishing time delay is recovered. The 
contributions from pairs of orbits with longer excursion 
times are much less relevant due to wave-packet spread- 
ing. 

Structured spectra, however, are only obtained if sev- 
eral pairs of orbits are taken into account. One should 
note that, if one employs orbits starting from (3,4) and 
vanishing time delays between the attosecond-pulse train 
and the IR field, there seems to be an improvement in 
the resolution of the spectra (c.f. blue line in Fig. El (a)). 
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FIG. 6: (Color online) Contributions of isolated pairs of 
orbits to the high- harmonic spectra for the same atomic 
and external- field parameters in Figs. 4 and 5, and delays 
(j> = 0, (j) = 0.257r, (j) — O.Stt and <j!> = 0.757r between the 
attosecond-pulse train and the low-frequency driving wave 
(Figs, inia), |Hl(b), |Hl(c) and|Sl(d), respectively). In the fig- 
ure, we only show the harmonic orders for which Q > Ip. The 
contributions from specific pairs of orbits are depicted in the 
same colors as those employed in Fig. 5, for depicting each 
corresponding pair. For comparison, the spectra computed 
with the ten shortest orbits and with orbit 3 to 10 are given 
as the thick lines in the figure. The numbers on the upper 
parts of the figure give the approximate kinetic energy of the 
electron upon return, in units of the ponderomotive energy. 



In |g, resolved harmonics have been related to select- 
ing a particular trajectory for the returning electron. In 
our case, however, it appears to be an interference effect 
between several orbits for the returning electron. The 
fast oscillations observed here are possibly not harmon- 
ics (due to the periodicity of the process), but are caused 
by the interference of different quantum orbits originat- 
ing from (launched from) the same half cycle. 

With regard to the above discussion, we should point 
out that the absence of discrete harmonics in the spec- 
tra displayed in this work is due to the fact that we are 
considering orbits starting only in the first cycle of the 
laser field, instead of over many cycles. Therefore, we do 
not have enough periodicity to obtain discrete harmonic 
peaks. This is however not a problem in our context, 
since we concentrate our investigations on how quantum- 
interference effects influence the main structure of the 
spectra. By considering a restricted range of electron 
start times, such effects can be more clearly seen and 
interpreted. 

The above-discussed behavior is very similar to that 
observed in 0. Figs. 14151 however, have been plotted 
for a set of time delays different from those employed 
in such a reference. Such discrepancies may be related 



to the fact that, in |g, only a finite and relatively small 
group of harmonics has been used, whereas, in our com- 
putations, we have employed the attosecond-pulse limit 
for an infinite set of harmonics. 

In fact, if the attosecond-pulse train is composed by a 
relatively small number of harmonics (such as in 6]), it is 
possible to define a main frequency for this train. Hence, 
it is reasonable to assume that electron is being ejected 
in the continuum by a photon of such a frequency, and 
thus with nonvanishing velocity. This will cause a split- 
ting of the trajectories into the "downhill" and "uphill" 
transitions, which correspond to the ejection of the elec- 
tron in- or opposite to- the direction of the laser field, 
respectively, and which will be selected according to the 
time delay. 

On the other hand, if the attosecond pulse train is 
composed by an infinite number of frequencies, the un- 
certainty relation will no longer allow one to define a main 
frequency for it. Consequently, the physical picture of an 
electron being released by a photon of a specific frequency 
and reaching the continuum with a constant velocity is 
no longer applicable. However, it is still possible to de- 
fine a specific time t' for which the electron is ejected, 
and determine the corresponding return times. 



IV. CONCLUSIONS 

In conclusion, we have shown that an attosecond-pulse 
train superposed to a strong, near-infrared laser field can 
be used for controlling above-threshold ionization (ATI) 
and high-order harmonic generation (HHG). This control 
is exerted by using the attosecond pulses, which are time 
delayed with respect to the infrared field, to constrain 
the instant in which an electronic wave packet is ejected 
in the continuum. Since this wave packet will either re- 
combine with its parent ion, generating high harmonics, 
or contribute to the above-threshold ionization spectrum, 
both phenomena are strongly affected by such pulses. 

Similar investigations, which led to similar effects to 
those in this paper, have been performed in em- 
ploying the numerical solution of the Schrodinger equa- 
tion. In this work, however, we go beyond such studies, 
and trace all features observed to quantum-interference 
effects between the possible pathways of the electronic 
wave packet. In particular, the model used in this paper, 
based on the Strong-Field Approximation (SFA) allows 
one to investigate the different ingredients involved in the 
computation of the transition amplitude, and how they 
are infiuenced by the attosecond-pulse train, in a detailed 
and transparent way. 

For ATI, we restrict ourselves to the direct electrons, 
i.e., those that reach the detector without rescattering 
with its parent ion. We consider an attosecond-pulse 
train composed by an infinite group of harmonics and we 
sum over all periods of the driving field. In this case, we 
observed an extension in the photoelectron energy if the 
attosecond-pulse train is in phase with the IR field (i.e.. 
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for delays (f) = nn). 

In j8(|, this has been attributed to the fact that, in this 
case, the electrons are ejected at a zero crossing of the IR 
field. Hence, they gain the maximal possible momentum 
from the field, namely p = ±2y/lJp. In this paper, we 
have shown that this mechanism does play a role, which, 
in our framework, is implicit in the non-trivial part of the 
action. However, in order to obtain the above-mentioned 
effect, the dipole matrix elements from a Is hydrogenic 
state are also a necessary ingredient. Physically, such 
matrix elements are related to the initial state where the 
electron is bound and, in the SFA, contain all the infor- 
mation about the atomic bound potential. Furthermore, 
we have shown that such an effect only occurs if the pho- 
toelectron energy is sufficiently far from the ionization 
threshold. In fact, for low ATI orders, maximal intensi- 
ties are obtained for cj) = {2n + 1)tt/2. This is also due 
to an interplay between the phase the electron acquires 
when it reaches the continuum and its initial bound state. 

Apart from such studies, we have observed that the 
ATI angular distributions are strongly dependent on this 
delay, and on the photoelectron energy. Near the ion- 
ization threshold, such distributions exhibit very pro- 
nounced maxima for backscattered electrons (alignment 
angle 6 = rnr) which leave at peak-field times (0 = 
(2n -I- l)7r/2). As the photoelectron energies increase, 
such maxima move towards the zero crossings of the elec- 
tric field {(p = 2mi). This behavior is in agreement with 
that observed for the ATI spectra. 

For HHG, depending on the time delay between the 
attosecond-pulse train and the IR field, we observed de- 
viations in almost 30% in the cutoff energy. Furthermore, 
the plateau structure is highly sensitive to the time delay 
(j). Indeed, we found that, for certain ranges of this pa- 
rameter, the low-plateau harmonics are one to two orders 
of magnitude stronger than those near its high-energy 
end. The precise sets of harmonics in such a double- 
plateau structure depend very much on cf). For our par- 
ticular parameters, the set of strong harmonics moved 
towards lower energies as 4> distanced itself from a zero 
crossing of the electric field. 

Our results can be traced back to the interplay of two 
sets of orbits for an electron returning to its parent ion. 
Specifically, the maximal return electron energy in the 
shortest set of orbits changes from to 1.8J7p, as the 
ejection time of the electron is varied. These orbits yield 
a group of very strong harmonics, whose extension can 
be manipulated by controlling the time delay <j). The re- 
maining harmonics are mainly determined by the second 
shortest pair of orbits, and are considerably weaker. 

In this context, a noteworthy aspect is that, in con- 
trast to ATI, for which the effects mentioned in this pa- 
per also depend on the initial electron state, in HHG they 
appear to be due to controlling electron orbits with an 
attosecond-pulse train. Indeed, whereas in ATI such ef- 
fects do not occur for initial localized bound states, i.e., 
in the broad gaussian approximation, for HHG they are 
already present in this case. 



The results for HHG also agree to a great extent with 
those reported in 6J. However, their interpretation, as 
well as the time delays for which they occur, are some- 
what different. In fact, in all features are explained by 
selecting a single orbit from the shortest pair, and by the 
splitting in this set of orbits due to the fact that the elec- 
tron is reaching the continuum with non- vanishing veloc- 
ity. Such differences are possibly due to the fact that, in 
16] , only a small group of harmonics has been employed 
for producing the attosecond-pulse train, while, in Eq. 

we consider the limit of an infinite number of har- 
monics. We therefore expect the agreement to increase if 
a larger number of harmonics is taken when constructing 
the attosecond-pulse-train. 

V. APPENDIX 

In this appendix, we explicitly provide the time- 
dependent expression adopted in Ref. 0|. In this case, 
the time-dependent dipole moment reads 

d{t) = i I dt' I d^fcd;(k + A(t))4(k + A(i')) 

E{t')eyiv[iS{t,t',\i)], (19) 
with the action 

5(t, t', k) = -i ^*[k + A(T)]2rfr - Ip{t - t'). (20) 

Alternatively to the procedure adopted in this paper, 
one may integrate Eq. I|19|) over the intermediate momen- 
tum using the saddle point equation (|15|l . This yields the 
equation 

h = '-^ J2i-ir h-Kt - 4 + (21) 

ri=0 

d:(k,t + A(t))4(ksi + A(4)) exp[i^(t, 4, k,0], 

with the action given by formula (20), to which a Fast 
Fourier transform routine is then applied. In Eq. (21), 
k,4 = ~{l/{t - t'J) dTA(r), and the sum extends 
over i'j < t. Due to the wave-packet-spreading effects 

(described by the {{2n/i{t ~t'^ + ie))^^^), the mam con- 
tributions to the sum in Eq. (21) come from few that 
are the closest to t. 
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